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The  research  performed  under  the  subject  grant  dealt  with  two  problems  in 
which  body  forces  lead  to  complex  flows  which  are  of  scientific  and 
technological  importance*  The  first  of  these  is  concerned  with  the  central 
unresolved  problem  of  flows  in  completely  confined  enclosures,  vie*,  that  it  is 
not  possible  to  predict  the  flow  pattern  a  priori  for  a  specified  geometric 
configuration  and  Imposed  boundary  conditions*  The  second  problem  is  an 
external  one  in  which  there  is  an  interaction  of  hydrodynamic  (Taylor-Gortler) 
and  thermal  (Raylelgh-B^nard)  instabilities*  The  details  of  the  work  done  are 
presented  below. 

Natural  Convection  in  Shallow  Enclosures 

Natural  convection  in  enclosures  is  encountered  in  many  technologies  such 
as  materials  processing,  crystal  growth,  solar  energy  receivers,  and  fluids 
storage  and  management*  The  associated  transport  phenomena  are  complex  and 
have  been  found  to  be  very  sensitive  to  the  geometric  configuration  and  the 
imposed  boundary  conditions*  As  a  result,  each  problem  has  to  be  solved 
separately*  Such  solutions,  however,  are  difficult  to  obtain  because  it  has 
not  been  possible  to  predict  a  priori  the  general  flow  pattern  for  a  given 
geometric  configuration  and  Imposed  boundary  conditions*  The  knowledge  of  the 
general  flow  pattern  is  not  only  Important  for  analytical  studies  but  also  for 
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numerical  and  experimental  investigations  because  it  indicates,  at  the  very 
least,  the  type  of  resolution  required  to  obtain  valid  results.  The  first  main 
objective  of  the  research  under  the  subject  grant,  therefore,  was  to  try  to 
resolve  this  central  problem  of  flow  pattern  prediction.  A  combined 
theoretical  and  experimental  program  was  carried  out  to  that  end.  Emphasis  was 
given  to  low  aspect  ratio  (shallow)  enclosures  because  of  recent  practical 
interest  in  such  geometries  and  because  all  the  appropriate  physics  common  to 
all  confined  natural  convection  is  contained  in  such  situations.  Thus,  the 
second  objective  was  to  obtain  detailed  data  for  shallow  enclosures.  A 
comprehensive  recent  survey  of  natural  convection  in  enclosures  la  given  by 
Ostrach  [1]. 

Experimental  Program: 

Experiments  in  shallow  enclosures  heated  at  the  ends  with  linear 

6 

horizontal  wall  temperatures  for  0.055  <  A  <  0.5,  27.7  <  Gr^  10  ,  and  0.72 
3 

Pr  <  1.38  x  10  are  described  by  Ostrach  et  al  [2],  where  A  is  the  height-to- 
length  (aspect)  ratio,  GrH  is  the  Grashof  number  based  on  the  height,  and  Pr  is 
the  Prandtl  number.  It  was  found  that: 

a.  The  overall  flow  patterns  are  predominantly  unicellular 

b.  The  flow  pattern  is  Inclined  relative  to  the  enclosure  geometry  for 
aspect  ratios  greater  than  approximately  0.1.  For  smaller  aspect 
ratios  it  becomes  essentially  parallel. 

c.  The  skewness  of  the  streamlines  and  flow  velocities  are  proportional 
to  the  Grashof  number  and  aspect  ratio. 

d.  Secondary  cells  are  observed  near  the  enclosure  ends  when  the  aspect 
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ratio  is  0.2  and  (Gr^PrA)  is  of  the  order  of  10  .A  representative 
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A  corresponding  comprehensive  study  of  the  heat  transfer  aspects  Is 
presented  by  Kaaotani  et  al  [3].  The  Influence  of  the  secondary  cells  os  the 
heat  transfer  Is  clearly  shown  therein,  and  the  inadequacy  of  analyses  that  do 
not  treat  the  flow  pattern  correctly  Is  demonstrated. 


Theoretical  Program: 

The  experimental  prograa  has  indicated  that  in  shallow  enclosures  the  flow 
pattern  can  be  parallel  or  skewed  and,  under  certain  conditions,  secondary 
cells  or  flow  subregions  can  occur.  Such  ambiguities  concerning  the  nature  of 
the  flows  have  been  found  for  other  configurations  and  are,  in  fact,  inherent 
to  all  internal  flows.  Since  there  was  no  way  to  predict  the  flow  pattern  with 
confidence  a  theoretical  program  was  initiated  to  develop  a  method  for 
obtaining  a  qualitative  picture  of  the  overall  flow  pattern  for  a  given 
geometric  configuration  and  imposed  boundary  conditions.  Because  of  the  more 
extensive  experimental  data  for  shallow  enclosures,  emphasis  was  given  to  that 
configuration.  The  first  phase  of  this  work  is  described  by  Lee  and  Ostrach 
[4].  The  essence  of  that  work  and  the  results  will  be  briefly  smmrited. 

In  low  aspect  ratio  rectangular  enclosures  the  flow  pattern  can,  on  a 
geometric  basis  alone,  be  essentially  divided  into  two  regiona,  vis.,  one,  the 
region  near  the  end  walls  which  is  referred  to  as  the  end  region  and  the  other, 
the  core  region,  is  the  region  outside  the  end  regions  and  it  is  bounded  by  the 
horizontal  boundaries. 

In  order  to  be  able  to  predict  the  core  flow  pattern  and,  thus,  the  entire 


flow  pattern  correctly  one  has  to  be  clear  as  to  what  physical  mechanisms 
pertain  to  each  region.  Only  in  this  way  can  an  appropriate  mathematical  model 
be  developed.  Since  the  flow  characteristics  can  be  different  and  coupled  in 

each  flow  region  it  is  not  possible  to  consider  all  the  important  physical  codes" 

id/or 
?peeial 


A 


□Q 


■echnliu  in  the  enclosure  from  the  basic  equations  which  apply  in  either  one 
region.  In  addition,  the  equations  that  describe  the  phenomena  of  Interest 
herein  are  non-linear  and  bidirectionally  coupled.  Any  s lap lif lea t ions  must, 
therefore,  be  carefully  asde  and  ad  hoc  assumptions  are  dangerous.  In  order  to 
ensure  that  an  assumption  made  in  one  equation  is  consistently  transmitted  to 
the  other  equations  a  somewhat  formal  procedure  is  employed  based  on  the  method 
of  multiple  scales.  Multiple  scales  are  introduced  to  glue  mathematical 
degrees  of  freedom  which  enable  physical  statements  to  be  made  properly 
concerning  the  Important  physical  mechanisms  in  the  enclosure.  Analytically 
the  method  of  multiple  scales  can  be  used  to  obtain  single  uniformly  valid 
expansions,  in  contrast  to  the  method  of  matched  asymptotic  expansions  which 
yields  two  expansions  that  aust  be  matched.  The  method  of  multiple  scales 
formally  yields  the  equations  valid  in  each  region  as  well  as  in  the 
interaction  region. 

With  the  available  degrees  of  freedom  obtained  from  the  multiple  scaling 
proper  physical  statements  (force  or  energy  balances)  can  be  nade  in  the 
equations  that  have  been  properly  normalized.  The  order  of  magnitude  of  each 
term  in  the  equaiton  can  then  be  estimated  from  the  dimensionless  parameters 
that  appear  as  coefficients  of  the  various  terms.  It  is  easy  to  extract  the 
equations  that  describe  the  core  flow  characteristics  from  the  general 
equations. 

The  core  flow  patterns  are  first  studied  in  a  global  sense  under  the 
implicit  assumption  that  there  exist  no  flow  subregions,  such  as  secondary 
cells.  On  that  basis,  the  geometric  length  scales  are  the  appropriate  ones  for 
the  core  flow  structure.  With  the  use  of  the  results  of  the  global 
considerations  detailed  analysis  is  then  made  to  determine  the  possibility  of 


f 


-5- 

f  low  subregions.  In  the  detailed  analysis  arbitrary  length  scales  which  are 

determined  from  the  physical  statements  replace  the  geometric  length  scales. 

The  proper  scaling  of  all  the  variables  are  obtained  in  this  way.  The  results 

to  be  presented  pertain  essentially  to  enclosures  with  adiabatic  horizontal 

surfaces  and  with  the  (vertical)  ends  at  different  temperatures.  For  some 

situations  other  horizontal  thermal  conditions  were  examined. 

The  results  are  summarized  in  Table  1.  The  basic  force  and  energy 

balances  in  each  region  are  indicated  thereon  as  are  also  the  proper  scaling 

for  the  characteristic  stream  function,  ^  r ,  and  the  horizontal  length  scale  of 

end  region,  6 x,  and  the  physical  conditions  for  which  the  various  core 

configurations  will  be  obtained.  Note  that  three  specific  situations  are 

analyzed  depending  on  where  the  buoyancy  (driving)  force  is  operative.  When 
2  2 

RajjA  <  1  and  PrRa^A  «  1  the  core  temperature  distribution  is  linear  and  the 
flow  is  driven  by  the  buoyancy  in  the  core.  Depending  on  the  magnitude  of  the 

Grashof  number  the  core  flow  will  be  parallel  when  there  is  no  flow  boundary 

2  2 
layer  (GrA  <  1)  and  nonparallel  when  there  is  a  flow  boundary  layer  (GrA  > 

1).  No  flow  subregions  are  expected  under  these  conditions. 

When  there  are  thermal  boundary  layers  near  the  ends  of  the  enclosure 
2 

(RagA  »  1)  the  flow  will  be  driven  by  the  buoyancy  acting  there.  Then 
depending  on  the  magnitude  of  Pr  the  core  flow  will  be  stagnant  in  the  mid- 
core,  distinct  horizontal  thermal  boundary  layers  will  exist,  and  the  core 
temperature  distribution  will  be  stratified.  Secondary  cells  can  be  expected 
for  large  Pr. 

In  the  intermediate  flow  regime  the  core  temperature  distribution  varies 
with  both  space  variables  so  that  the  buoyancy  force  acts  throughout  the 
enclosure  and  there  is  an  interaction  between  the  core  and  end  flows.  In  this 
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cu«  both  parallel  and  non-parallel  core  flows  can  occur  under  the  conditions 
specified  and  nultlcell  patterns  seen  possible. 

Conparison  of  the  predictions  of  the  core  flow  patterns  show  good 
agreement  with  existing  experimental  data  for  Pr  1.  For  Pr  <  1  there  is 
insufficient  data  for  a  conparable  configuration.  The  global  prediction  of  the 
core  flow  patterns  is  generally  satisfactory.  Further  work  however,  will  be 
required  to  define  the  flow  subregions  better.  The  method  of  analysis  has 
clarified  much  of  the  physics  of  such  problems  and  can  be  generalized  to  other 
configurations.  It  thus  seems  like  a  good  first  step  in  resolving  the  central 
problen  that  has  confronted  workers  on  this  subject  since  its  inception. 

Combined  Thermal  and  Hydrodynamic  Instability  in  the 
Boundary  Layer  along  a  Curved  Plate 

This  work  is  described  in  detail  by  Lin  et  at  [5]  so  that  only  a  summary 
of  the  work  is  presented  herein. 

The  Important  dimensionless  parameters  of  the  problem  are  G  (GOrtler 
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number),  Gr  (Grashof  number)  and  Pr  (Prandtl  number).  The  ratio  Gr/G  is 

important  in  the  present  study  because  it  represents  the  relative  importance  of 

buoyancy  to  centrifugal  forces  on  the  behavior  of  Gortler  vortices.  The  ranges 

of  the  parameters  in  the  present  experiments  are  Pr  ■  .7  (air),  G  <  8.1  and  0  < 
2 

Gr/G  £  8.5.  The  velocity  and  temperature  distributions  in  the  boundary  layers 
are  measured  in  detail  by  a  hot-wire  and  a  thermocouple  probe.  Based  on  the 
data  the  effects  of  destabilising  heating  on  Gortler  vortices  are  delineated. 
The  data  takan  under  Isothermal  conditions  (Gr  ■  0)  and  some  data  with  heating 
were  given  in  the  previous  reports.  As  reported,  Gortler  vortices  causm 
spanwise  variations  of  the  velocity  and  temperature  distributions. 


An  iaportnnt  quantity  to  characterise  the  strength  of  Gortler  vortices  is 
the  amplitude  of  the  spanvlse  velocity  or  temperature  variations*  Under 
isothermal  conditions  it  is  found  that  the  amplitude  of  the  spanvlse  velocity 
variation  is  closely  approximated  by  the  relation 

gj  ■  exp  [•  (G  -  Gc)] 

where  Gc  is  the  critical  Gortler  number,  Ao  the  amplitude  at  G  ■  Gc,  and  a  • 
.41.  The  above  correlation  curve  is  compared  with  the  experimental  data  in 
Fig.  2.  With  heating,  the  amplitude  of  the  spanvlse  temperature  variation  is 
found  to  be  expressed  as 

g£  -  exp  [a  (C  -Gc)] 

2  1/2 

where  C  is  the  combined  parameter  C  -  (G  +  f  Gr)  ,  and  f  represents  the 
relative  importance  of  buoyancy  to  centrifugal  forces,  f  is  found  to  depend  on 
the  free  stream  velocity.  The  correlated  curves  are  compared  with  the 
experimental  data  in  Fig.  3.  Comparison  of  the  above  two  expressions  shows 
that  one  important  effect  of  heating  is  to  Increase  the  strength  of  the 
vortices. 

Another  Important  effect  of  heating  is  that  it  enhances  the  non-linear 
development  process  of  the  vortices.  Fig.  4  shows  the  effect  of  heating  on  the 
spanvlse  temperature  distributions.  As  seen  in  the  figure,  with  increasing 
heating  the  temperature  distribution  in  the  vara  region  becomes  spiky,  which 
means  that  the  upwash  region  of  the  vortices  tends  to  be  concentrated  in  a 
narrow  region,  but  the  temperature  distribution  in  the  cold  region  suggests 
that  the  dovnvash  region  of  the  vortlce  tends  to  be  more  spread  out.  This  is 
the  non-linear  effect  Investigated  in  the  past  under  isothermal  conditions. 


The  present  study  shows  that  the  non-linear  effect  is  enhanced  hy  heating. 
Also  due  to  the  enhancement  the  vortices  start  to  meander  in  the  spanwlse 
direction,  which  causes  temperature  fluctuation  with  time.  The  temperature 
fluctuating  levels  at  three  spanwlse  sections  are  shown  in  Fig.  5.  Sections  V 
and  P  correspond  to  the  upwash  region  and  the  downwash  region  of  the  vortices, 
respectively.  Section  C  is  in  the  middle  of  the  two  sections.  Since  the 
temperature  fluctuation  is  mainly  caused  by  an  unsteady  meandering  motion  of 
the  vortices,  the  fluctuating  level  is  smallest  in  the  downwash  region  because 
of  relatively  uniform  temperature  in  the  region  as  described  above,  and  largest 
in  the  upwash  region  due  to  relatively  non-uniform  spanwlse  temperature 
variation  in  the  region. 

The  conclusions  drawn  from  the  present  experimental  study  are  as  follows. 

1.  Within  the  ranges  of  the  par ante ter  Gr/G^  studied  herein  (Gr/G*  <  8.5)  the 
wavelength  of  the  vortice  remains  unchanged.  The  strength  of  the  vortices  is 
enhanced  by  heating. 

2.  The  amplitude  of  the  vortices  Increases  almost  exponentially  with  the 
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combined  parameter  (G  +  f  Gr)  where  f  is  found  to  be  .3  <  f  <  .4,  until  the 
non-linear  effects  become  important. 

3.  Heating  enhances  the  non-linear  effect.  The  warm  upwash  region  of  the 
vortices  tends  to  be  concentrated  in  a  narrow  region  and  the  cold  downwash 
region  tends  to  spread  out  as  Gr  increases  while  G  is  fixed,  amd  as  Gr  is 
increased  further,  the  vortex  rolls  start  to  meander,  which  causes  temperature 
fluctuations  with  time  at  a  fixed  point. 

4.  The  non-linear  effects  -  the  vortex  shape  distortion  and  the  meandering  of 

2 

the  vortex  rolls  -  become  noticeable  at  a  smaller  value  of  the  parameter  G  +  f 

2 

Gr  as  the  value  of  the  parameter  Gr/G  is  increased. 
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